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Abstract 



We present a detailed and systematic analysis of the coarse-grained, nonequi- 
librium dynamics of a scalar inflaton field coupled to a fermion field in the 
latter stages of the reheating period of inflationary cosmology. We derive 
coupled nonperturbative equations for the inflaton mean field and variance at 
two loops in a general curved spacetime, and show that the equations of mo- 
tion are real and causal, and that the gap equation for the two-point function 
is dissipative due to fermion particle production. We then specialize to the 
case of Minkowski space and small- amplitude inflaton oscillations, and derive 
the perturbative one-loop dissipation and noise kernels to fourth order in the 
Yukawa coupling constant; the normal-threshold dissipation and noise kernels 
are shown to satisfy a zero-temperature fluctuation-dissipation relation. We 
derive a Langevin equation for the dynamics of the inflaton zero mode, and 
show that its variance is non-negligible during reheating. Stochastic evolution 
of the inflaton may have important astrophysical consequences. 
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I. INTRODUCTION 



In an earlier paper [|I| we studied the effect of parametric particle creation on the infla- 
ton dynamics in the post-inflation, "preheating" stage. It was shown that analysis of the 
dominant physical processes during the early stages of reheating necessitates consideration 
of self-consistent back reaction of the infiaton field variance on the inflaton mean field and 
quantum modes. This is a nonperturbative effect, and its description requires a consistent 
truncation of the Schwinger- Dyson equations Assuming initial conditions conducive 

to efficient parametric particle creation, the end state of the regime of parametric particle 
creation consists of a large infiaton variance (i.e., on the order of the tree level terms in 
the infiaton's effective mass). Thus, both the infiaton mean field and variance should be 
treated on an equal footing. This requires a two-particle-irreducible (2PI) formulation of 
the effective action which is a subclass of the master effective action 0. 

During the later stages of reheating, the dynamics of the infiaton field is thought, in the 
case of unbroken symmetry, to be dominated by damping from fermionic particle creation 
0]. This stage of infiaton dynamics is the subject of this paper. We consider a model 
consisting of a scalar infiaton field (with X(j)'^ self-coupling) coupled to a spinor field ip via 
a Yukawa interaction and we attempt to present as complete and rigorous a treatment as 
mandated by the self-consistency of the formalism and the actual solvability of the equations. 
Thus, we adopt a closed-time-path (CTP), two-particle-irreducible (2PI), coarse-grained 
effective action (CGEA) to derive the dynamics of the infiaton field. We have explained the 
significance of CTP and 2PI in our earlier papers p|-P,p|,p!0| and their relevance to the study 
of infiaton dynamics in Ref. 0]. Here, an added feature of an open system is introduced: we 
wish to include the averaged effect of an environment on the system, and a useful method 
is via the coarse-grained effective action |[TT|-[r5|. Let us explain the rationale for this. 

Coarse- Grained Effective Action: In inflationary cosmology at the onset of the reheating 
period, the infiaton field's zero mode typically has a large expectation value, whereas all 
other fields coupled to the infiaton, as well as infiaton modes with momenta greater than 



the Hubble constant, are to a good approximation in a vacuum state [^. This suggests 
imposing a physical coarse-graining in which one regards the infiaton field as the system, 
and the various quantum fields coupled to it as the environment. From the closed-time-path, 
coarse-grained effective action (CTP-CGEA) []TT]-0 derived in Sec. || below, one obtains 
effective dynamical equations for the infiaton field, taking into account its effect on the 
environment, and back reaction therefrom. For the present problem, the system consists 
of the infiaton mean field and variance, and the environment consists of the spinor field(s) 
coupled to the system via a Yukawa interaction. 

We wish to emphasize here a subtle yet important distinction between the system- 
environment division in nonequilibrium statistical mechanics and the system-bath division 
assumed in thermal field theory. In the latter, one assumes that the propagators for the bath 
degrees of freedom are fixed, finite-temperature equilibrium Green functions, whereas in the 
case of the CTP-CGEA, the environmental propagators are slaved (in the sense of [Q) to 
the dynamics of the system degrees of freedom, and are not fixed a priori to be equilibrium 
Green functions for all time. This distinction is important for discussions of fermion parti- 
cle production during reheating, because it is only when the infiaton mean field amplitude 
is small enough for the use of perturbation theory, that the system-bath split implicit in 
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thermal field theory can be used. Otherwise, one must take into account the effect of the 
inflaton mean field on the bath (spinor) propagators. 

Earlier work: The first studies of particle production during reheating in inflationary 
cosmology were where reasonable estimates of particle production were made, but 



back reaction effects were not addressed. The earliest studies of fermionic particle creation 
during reheating used time-dependent perturbation theory to compute the imaginary part 
of the self-energy for the zero mode of the inflaton field, which was related to the damping 
parameter in a friction-type phenomenological term in the equation of motion for the inflaton 
zero mode p!B| , pT| , pi?| , P^p3| . In these studies, it was assumed that the effect of fermionic 
particle production on the djTiamics of the inflaton zero mode is that of a Tcf) friction term. 
However, it has been shown |25|,[7| j26|j27|j2^j28|j2l|j6^ that this assumption is not tenable for a 
wide variety of field-theoretic interactions and initial conditions. Rather, the effect of back 
reaction from particle production must be accounted for by deriving the effective evolution 
equation for the quantum expectation value of the inflaton zero mode, where the dynamics 
of the degrees of freedom of the produced particles are either coarse-grained (as in Sec. |T| 
of this paper), or accounted for through self-consistent coupled equations (as in Sec. y). 
In general, particle creation leads to a nonlocal dissipation term in the inflaton mean field 
equation, and it is only under rather idealized conditions and specialized cases that one 
can expect the dissipation kernel to approach a delta function (i.e., a friction term) p9-32 



Therefore, ||16| , |22| , |23| missed the time-nonlocal nature of fermion particle production and its 
effect on the dynamics of the inflaton field. In addition, these studies computed the self- 
energy in flat space, thereby neglecting the effect of curved spacetime on fermion production, 
and did not examine the stochastic noise arising from the coarse-graining of the fermion field. 

In addition, most early studies of fermion production during reheating, in using time- 
dependent perturbation theory to compute the vacuum particle production rate |p!B| , P^ , p5| , 
did not include the effect of back reaction of the produced fermion particles on the particle 
production process itself. In |^T|, the effect of a thermal initial fermion distribution on 
the particle production process was investigated (and a Pauli blocking effect was shown), 
but their analytic derivation of the Pauli blocking effect involves the same perturbative 
expansion (i.e., system-bath split) described above, and therefore does not incorporate the 
effect of the produced fermion quanta on the particle production process. In order to take 
this effect into account, it is necessary to include the effect of the time- varying inflaton 
mean field in the equation of motion for the spinor propagator, which amounts to a coarse- 
graining of the fermion field, in the system-environment sense, as described above. The 
perturbative amplitude expansion of the effective inflaton dynamics, in contrast, amounts 
to a system-bath coarse-graining which does not include this back reaction effect. 

More recent studies of fermion production during reheating |]33| , |3^j2T[| obtained dynamical 
equations for the inflaton mean field using a one-loop factorization of the Lagrangian, and 
solved them numerically. However, these studies were carried out in flat space, and because 
they studied only the dynamics of the inflaton mean field (and at one loop), their analysis 
did not take into account the back reaction of the inflaton variance on the fermion quantum 
modes, nor the back reaction of particle production on the dynamics of the inflaton two-point 
function. The importance of the curved spacetime effect has been addressed in a previous 
paper , and we will discuss below the importance of treating the inflaton quantum variance 
on equal footing with the mean field. 
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Present work: In this study, we wish to describe the late stages of the reheating period, in 
which the damping of the inflaton mean field is dominated by fermionic particle production; 
our starting point is the end of the preheating period (in which the inflaton dynamics was 
dominated by back reaction from parametric particle creation, as discussed in our earlier 
paper Because the inflaton variance (y?^) can (for sufficiently strong self-coupling) be on 
the order of the square of the amplitude of mean field oscillations at the end of preheating 
(in the case of unbroken symmetry) it is necessary to treat the inflaton mean field and 
variance on an equal footing in a study of the subsequent effective dynamics of the inflaton 
field. This requires a two-loop, two-particle-irreducible formulation of the coarse-grained 
effective action. At two loops, both the inflaton mean field and the inflaton variance couple 
to the spinor degrees of freedom, and are damped by back reaction from fermion particle 
production; all the previous studies mentioned above, in using the IPI effective action, 
missed this possibly important effect. 

In addition to having a large variance, the inflaton amplitude at the end of the pre- 
heating period may be large enough that the usual perturbative expansion in powers of the 
Yukawa coupling constant is not valid [see Eq. ( p.l| ) below], in which case a nonperturba- 
tive derivation of the inflaton dynamics is required. In the construction of the CTP-2PI, 
coarse-grained effective action below, the spinor propagators obey one-loop dynamical equa- 
tions in which the inflaton mean field acts as an external source. Studies of fermion particle 
production during reheating which relied on the use of perturbation theory in the Yukawa 



coupling constant |16,r8|JT7|JT9|j24--p4| therefore do not apply to the case of fermion particle 



production at the end of preheating with unbroken symmetry, when the Yukawa coupling 
is sufficiently large. The dynamical equations derived in Sec. |II| below for the inflaton mean 
field and variance are applicable even when, as may be the case, the inflaton mean field 
amplitude is large enough that a perturbative expansion in powers of the Yukawa coupling 
is not justified. 

Although, as discussed above, a proper treatment of the early stage of fermion produc- 
tion during reheating, starting at the end of preheating with unbroken symmetry, should 
in principle employ the CTP-2PI-CGEA to obtain coupled equations of motion for the in- 
flaton mean field and variance, at very late times the inflaton mean field and variance will 
have been damped sufficiently (due to the dissipative mechanisms derived below in Sec. |lD 
that the perturbative IPI effective action will yield a qualitatively correct description of 
the inflaton mean-field dynamics. While curved spacetime effects should in principle be 
incorporated self-consistently for a quantitative calculation of the reheating temperature in 
a particular inflationary model (as discussed in |jl[), for a general discussion of dissipative 
effective dynamics of the inflaton mean field in the case of weak cosmic expansion (where 
the Hubble constant H is much smaller than the frequency of inflaton oscillations), it is 
reasonable to neglect curved spacetime effects in computing the spinor propagators. There- 
fore in Sec. [TTTI , we derive the perturbative, flat-space CTP-IPI-CGEA to fourth order in 
the Yukawa coupling constant, and obtain an evolution equation for the inflaton mean field 
with nonlocal dissipation. 

Another new feature of our work obtainable only from the stochastic approach adopted 
here is the derivation, in Sec. 0, of a Langevin equation for the inflaton mean field, with 
clear identification of the dissipation and noise kernels from the CGEA. We have calculated 
the energy dissipated and the fluctuations in the energy. From the latter we obtain the range 
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of parameters where the conventional "mean-field" approach breaks down. We believe the 
methodology presented here provides a better theoretical framework for the investigation 
of phase transitions in the early universe, as exemplified by our treatment of reheating in 
infiationary cosmology. 

Organization and Notation: Our work is organized as follows. In Sec. |I|, we derive the 
coupled equations of motion for the infiaton mean field and variance, in a general curved 
spacetime, including diagrams in the coarse-grained CTP-2PI effective action of up to two- 
loop order. In Sec. |T|, we specialize to Minkowski space and small mean-field amplitude, 
and obtain a perturbative mean field equation including dissipative effects up to fourth 
order in the Yukawa coupling constant. In Sec. we examine the dissipation and noise 



kernels obtained in Sec. |lTl|, and show that they obey a fluctuation-dissipation relation. We 
then derive a Langevin equation which self-consistently includes the effect of noise on the 
dynamics of the inflaton zero mode. We summarize our results in Sec. 0. 

Throughout this paper we use units in which c = 1. Planck's constant h is shown explic- 
itly (i.e., not set equal to 1) except in those sections where noted. In these units, Newton's 
constant is G = ^Mp~^, where Mp is the Planck mass. We assume a four-dimensional, 
globally hyperbolic, orientable spacetime manifold, and follow the sign conventions^ of Bir- 
rell and Davies for the metric tensor gf^,^, the Riemann curvature tensor Rfj,urjpi and the 
Einstein tensor G^^. We use greek letters to denote spacetime indices. The beginning latin 



letters a, 6, c, rf, e, / indicate the time branch (see Ref. |T^, Sec. II). The Einstein summation 
convention over repeated indices is employed. 

II. COARSE-GRAINED INFLATON DYNAMICS IN CURVED SPACETIME 

In this section, we present a flrst-principles derivation of the nonequilibrium, nonpertur- 
bative, effective dynamics of a scalar inflaton fleld coupled to a spinor fleld via a Yukawa 
interaction, in a general, curved classical background spacetime. The use of the Schwinger- 
Keldysh closed-time-path (OTP) formalism allows formulation of the nonequilibrium dy- 
namics of the inflaton from an appropriately deflned initial quantum state. The evolution 
equations for the inflaton mean fleld and variance are derived from the two-loop, closed-time- 
path (CTP), two-particle- irreducible (2PI), coarse-grained effective action (CGEA). As the 
name suggests, there are two approximations of a statistical mechanical nature. One is the 
coarse graining of the environment — here the inflaton fleld is the system and the fermion 
fleld is the environment []TT|. The other refers to a truncation of the correlation hierarchy 



for the inflaton fleld — the two-particle-irreducible effective action. This formulation re- 
tains the inflaton mean fleld and variance as coupled dynamical degrees of freedom. Back 
Reaction of the scalar and spinor fleld dynamics on the spacetime is incorporated using the 
semiclassical Einstein equation, which follows from functional differentiation of the effective 
action with respect to the metric. It is shown that these dynamical equations are both real 



^In the classification scheme of Misner, Thorne and Wheeler [Q, the sign convention of Birrell 



and Davies |36| is classified as (+,+,+) 
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and causal, and the equations are cast in a form suitable for implementation of an explicit 
curved-spacetime renormalization procedure. 



A. The model 



We study a model consisting of a scalar field (the inflaton field) which is Yukawa- 
coupled to a spinor field ip, in a curved, dynamical, classical background spacetime. The 
total action 

S[<j>,^,^,gn = SV]+S'[<l>,^,^,9n (2.1) 
consists of a part depicting classical gravity, S'^[g'^^], and a part for the matter fields, 

S^[<P,ij,ij,gn = S'^[<p,gn + S^[iJ,iJ,gn + S''[<P,i^,i^,9n (2-2) 
whose scalar (inflaton), spinor (fermion), and Yukawa- interaction parts are given by 



S^l^^tfj^g'^''] 



d^X^/^ 



(V-T^V^V - (V^^)7^^) - 1^^^ 



(2.3a) 
(2.3b) 
(2.3c) 



For this theory to be renormalizable in semiclassical gravity, the bare gravity action S'^[g^'^] 
of Eq. ( |2.1| ) should have the form 



S''[(P,^P,ij,gn = -f / d^x^cf^il:. 



SV] = TTT^ I d'x^ [R - 2A, + cR^ + hR^^R^^ + aR^^^' R^p,s\ 
ibvrG J 



(2.4) 



In Eqs. ( p.3a| )- (|2.3c|) , m is the scalar field "mass" (with dimensions of inverse length); ^ is 
the dimensionless coupling to gravity; is the spinor field "mass," with dimensions of inverse 
length; □ is the Laplace- Beltrami operator in the curved background spacetime with metric 
tensor g^^] is the covariant derivative compatible with the metric; is the square root 
of the absolute value of the determinant of the metric; A is the self-coupling of the inflaton 
field, with dimensions of 1/v^; and / is the Yukawa coupling constant, which has dimensions 
of l/\/h. In Eq. (|2.4| ), G is Newton's constant (with dimensions of length divided by mass); 
R is the scalar curvature; R^y is the Ricci tensor; Ra/S'ys is the Riemann tensor; a, b, and c 
are constants with dimensions of length squared; and Ac is the cosmological constant, which 
has dimensions of inverse length-squared. The curved spacetime Dirac matrices 7^ satisfy 
the anticommutation relation 



{7^7n+ = 2^?'^ap, 



(2.5) 



in terms of the contravariant metric tensor g^'^. The symbol l^p denotes the identity element 
in the Dirac algebra. 



6 



Due to the generalized Gauss-Bonnet theorem the constants a, b, and c are not all 



independent in four spacetime dimensions; let us therefore set a = 0. It is assumed that 
there is a definite separation of time scales between the stage of "preheating" discussed 
in our earlier work and fermionic particle production, which is our primary focus in 
this work. However, this does not imply that perturbation theory in the Yukawa coupling 
constant / is necessarily valid, which would require that condition (|3.1| ) (defined in Sec. |T| 
below) be satisfied. In addition, the fermion field mass /i is assumed to be much lighter than 
the inflaton field mass m, i.e., the renormalized parameters m and fi satisfy m ^ fi. 

B. Closed-time-path, coarse-grained effective action 

We denote the quantum Heisenberg field operators of the scalar field and the spinor 
field ip by $h and \E'h5 respectively, and the quantum state^by For consistency with the 
truncation of the correlation hierarchy at second order, we assume $h to have a Gaussian 
moment expansion in the position basis ||39|, in which case the relevant observables are the 



scalar mean field 

4>{x) = {n\^^{x)\n), (2.6) 

and the mean-squared fluctuations, or variance, of the scalar field 

{n\<^l{x)\n) - (^]|$H(x)|^])2 = {Q\y,,\x)\n), (2.7) 

where the last equality follows from the definition of the scalar fluctuation field 

¥Ph(x) = $h(x) -0(x). (2.8) 

As discussed above, at the end of the preheating period, the inflaton variance can be as large 
as the square of the amplitude of mean-field oscillations. On the basis of our assumption 
of separation of time scales in Sec. |I1 A| , and the conditions which prevail at the onset of 
reheating, the initial quantum state \ fl) is assumed to be an appropriately defined vacuum 
state for the spinor field. 

The construction of the CTP-2PI-CGEA for the (pipip theory in a general, curved, back- 
ground spacetime closely parallels the construction of the CTP-2PI effective action for the 
0{N) model in curved spacetime discussed in an earlier paper |T0|. Within the spacetime 
manifold (whose dynamics must be determined self-consistently through the semiclassical 
gravitational field equation), let M be defined as the past domain of dependence of a Cauchy 
hypersurface T,^,, i.e., M = D_(S^), where has been chosen to be far to the future of any 
dynamics we wish to study. We now define a "OTP" manifold Ai as the quotient space 

^^(Mx {+,-})/-, (2.9) 



^Although in this case the particular initial conditions constitute a pure state, this formalism can 
encompass general mixed-state initial conditions 0. 
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in terms of the discrete set {+, — } which labels the time branch. The volume form sm on 
A4 is defined to be em on the plus (+) branch and —sm on the minus (— ) branch. The 
restrictions of a function (p, defined on A^, to the + and — time branches are denoted by 
and 0_, respectively. We can then define a matter field action on A4, 

5^ [0_ , ^_ , ^_ , g^^' ; 0+ , ^+ , ^+ , ] = [0+ , ^+ , , ^7r] - 5^ , ^- , ^- , gn , (2. lo) 

where the spacetime integrations in 5*^ are now over M only. We use the symbol to 
distinguish it from the action 5*^ on M. Let us also simplify notation by suppressing time 
branch indices in the argument of functionals on A4, i.e., 

S'[^,^,^,gn^s^[<P-,^-,i^-,9^'';<P+,i^+,i^+,9n (2-11) 

Let us also define the functional on M. by 

5-[0,^,^,^7n = ^"[0+,^+,^+,^?+1-^"[</>-,^-,^-,^7'^1, (2-12) 

in analogy with Eq. ( p.lO|) . For a function on the restrictions of to the + and — 
time branches are subject to the boundary condition 



(2.13) 



at the hypersurface S*. The gravity action 5*°, promoted to a functional on A4, takes the 
form 

S^[gt\gn=S'^[g'r]-S^[gn, (2.14) 

where the range of spacetime integration in on the right-hand side of Eq. (|2.14|) is 
understood to be over M. 

To formulate the CTP-2PI-CGEA, our first step is to define a generating functional for 
n-point functions of the scalar field, in terms of the initial quantum state \Q) which evolves 
under the influence of a local source J, and a non-local source K coupled to the scalar field 
(in the interaction picture with the external sources being treated as the "interaction"). This 
generating functional depends on both J and K, as well as the classical background metric 
g'^'^. In the path integral representation, the generating functional Z[J, K, g^^] takes the form 
of a sum over scalar field configurations and complex Grassmann-valued configurations ■0 
on the manifold 

Z[J,K,g>^''] = j^^ D^_Di,,Dij_D^^Dij+Dij+exp ^ (^^[0, ^, ^, ^^'^j 

' 4^ / ~„ab 



M 



A f [ rfVv^c'^V'^i^ac(a:,x')0b(x)0,(x'; 



(2.15) 



where Ja{x) is a local c- number source and Kab{ nonlocal c-number source. The 

subscript CTP on the functional integral denotes a summation over field configurations 0±, 
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il!±, and which satisfy the boundary condition (|2.13|) . The latin indices a, b, c, . . . , have 
the discrete index set {+,—}, and denote the time branch (^H- The boundary conditions 
on the functional integral of Eq. ( |2.15|) at the initial data surface determine the quantum 
state The CTP indices have been dropped from g^" for ease of notation; it will be clear 
how to reinstate them |T^ in the two-loop CTP-2PI effective action shown below in Sec. fT|. 
The two-index symbol c"'' is defined by 



1 if a = 6 = +, 
— 1 if a = 6 = — , 
otherwise. 

The generating functional for normalized n-point functions is 

W[J,K,g^''] = -th\nZ[J,K,g'^% 
in terms of which we can define the classical scalar field on 



6Jb{x) 

and the scalar two-point function on A4 in the presence of the sources Ja and K^b, 



(2.16) 



1 1 6W[J,Kg''' 



- <pa{x)(t)b{x') 



(2.17) 



(2.18) 



(2.19) 



where the JK subscripts indicate that and G are functionals of the Ja and Kab sources. 
In the limit Ja = Kab = 0, the classical field is the same on the two time branches, 



0^ 



J=K=0 



J=K=0 



(2.20) 



and it is equivalent to the mean field (p. In the same limit, Gab becomes the CTP propagator 
for the fluctuation field defined in Eq. (|2.8| ), 

(2.21a) 
(2.21b) 
(2.21c) 
(2.21d) 



hG++{x,x')\j^j^^Q = {Q\T{if^{x)if^{x'))\Q), 
hG-^{x,x')\j^j^^Q = {n\T{ipji{x)ipn{x'))\n), 

hG+^{x,x')\j^j^^Q = {n\(p^{x')ipj,{x)\n), 
hG^+{x',x)\j^j^^Q = {n\ip^{x)ip^{x')\n), 



where T denotes temporal ordering and T denotes antitemporal ordering. In the coincidence 
limit0 x' = X, all four components ( p.21a]) -( p.21d|) are equivalent to the variance {f-a^) 
defined in Eq. ( |2.7D . Provided we can invert Eqs. ( p^.lSj ) and ( |2.19| ) to obtain Ja and Kab in 
terms of (pa and Gab, the CTP-2PI effective action can be defined as the double Legendre 
transform (in both Ja and Kab) of W[J, K, g'^"], 



^The variance ((^(x)^) is divergent in four spacetime dimensions, and should be regularized using 
a covariant procedure [36,4^]. 
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T[<P,G,g^'']=W[J,K,g^'']- / d^x^gc^' Ja{x)<P,{x) 

J M 

-If d^x^ ! dVv^c'^V'='i^ac(x,x') \hGM{x,x')+^,{x)^a{x')\. {2.22) 

2 J M J M L ^ 

The inverses of Eqs. ( p. 18 ) and ( p.l9| ) can be obtained by functional differentiation of 
Eq. ( p.22|) with respect to 0^, 



1 m^,G,g' 



1 



-c'^Vb(x),, - -c'^V'^ / d^x'^[KUx,x'),, + Kd,{x\x),a]<P,{x'), 



M 



and with respect to Gah, 

1 1 5V[^,G,g^^-] 



h 



^ 6Gab{x,x') 



c^^d'^K^Ax^x' 



(2.23) 



(2.24) 



where the ipG subscript indicates that Kab and J a are functionals of (pa and Gab- Inserting 
Eqs. (|2.23|) and (|2.24| ) into Eq. ( p.22|) yields a functional integrodifferential equation for the 
CTP-2PI effective action in terms of and G only, so the JK subscripts can be dropped. 
It is useful to change the variable of functional integration to be the fluctuation field about 
0a, defined by 



■'a fa- 



(2.25) 



Performing the change-of-variables Dcf) —>■ Dip, the equation for F is 



mG,g 



flU] 



M 



d X d x' 



M 5Gba{,x\x) 



-Gab{x,x') 



-ih\ii{ I Dip+D'ilj+Dip+Dip_D'ilj_D'ilj_exp^-i^S^[ip + (f),'^,'ilj,g^"'] 



4 6r[<p,G,gn 

d X -ifa 



M 



-I d^X I d^X' 



M 



M 



jmG,g^^' 
6Gba{x',x) 



^-Mx)Mx')) \, (2.26) 



which has the formal solution 



r[0,G,(?n =5'^[0] - ylndetG'a5-^mndetF,fc + r2[0,G] 



+ '41 d'x^/^ [ d^x'^/^A''\x',x)Gab{x,x'), (2.27) 

where is the second functional derivative of the scalar part of the classical action S"^, 
evaluated at 0, 



iA'^Hx^x') 



()a(x)(50b(x') J 



X 



c-b^n, + + ^R{x)) + c''"''-(l),ix)Mx) 



6{x — x') 



(2.28) 
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The symbol Fab denotes the one-loop CTP spinor propagator, which is defined by 

Fai,ix,x')=B-,\x,x'), (2.29) 
where we are suppressing spinor indices, and the inverse spinor propagator B"''' is defined by 



iB^Hx^x) = ^ 



-9 



(2.30) 



[c'^'itYV^ - /i) - c»''70c(a;')) Six' - x)^l.,. (2.31) 



-9 

It is clear from Eq. ( |2.31|) that the use of the one-loop spinor propagators in the construction 
of the CTP-2PI-CGEA represents a nonperturbative resummation in the Yukawa coupling 
constant, which (as discussed above) goes beyond the standard time-dependent perturbation 
theory. The boundary conditions which define the inverses of Eqs. ( p.28| ) and ( p. 311) are the 
boundary conditions at the initial data surface in the functional integral in Eq. ( p. 151 ), which 
in turn, define the initial quantum state The one- loop spinor propagator Fab is related 
to the expectation values of the spinor Heisenberg field operators for a spinor field in the 
presence of the c-number background field 0, 

nF++(x,x')|^^=^_: = {n\Ti^M^,ix')m, (2.32a) 

nF__(x,x')|,^=,_=, = {n\fi^^ix)^^ix')m, (2.32b) 

h F^_{x,x%^^,_^, = -{n\^^ix')^M\^), (2.32c) 
hF_+{x,x%^^^_^^={n\^^{x)^^{x')\n), (2.32d) 

where the spinor Heisenberg field operators obey the equations 

(^7^V^-/W-/0)^ = O, (2.33a) 
(-z7'^V^-/i-/0)^ = O. (2.33b) 

The CTP spinor propagator components satisfy the relations (valid only when 0+ = 0_ = 0) 

F++{x,x'y = F^^{x',x), (2.34a) 

F__{x,xy = F++{x',x), (2.34b) 

F_+(x, x'y = F.+ {x\ x), (2.34c) 

F+_(x,x')^ = F+_(x',x), (2.34d) 

The functional r2[0, G] is defined as —ih times the sum of all vacuum diagrams drawn 
according to the following rules: 

1. Vertices carry spacetime {x G M) and time branch (a G {+,—}) labels. 

2. Scalar field lines denote hGab{x,x'). 

3. Spinor lines denote the one-loop CTP spinor propagator hFab{x, x') (spinor indices are 
suppressed), defined in Eq. ( ^^.291 ). 



11 



4. There are three interaction vertices, given by iS^/h, which is defined by 

S'[(j),(p,i>,ilj] = S'[(j)+,(p+,^|J+,^p+] - S"[0_,v2-, V"-,^-], (2.35a) 



(2.35b) 



where we have followed the notation of Eq. ( ^^.121 ). 



5. Only diagrams which are two-particle-irreducible with respect to cuts of scalar lines 
contribute to 

The distinction between the CTP-2PI, coarse-grained effective action which arises here, and 
the fully two-particle- irreducible effective action (2PI with respect to scalar and spinor cuts), 
is due to the fact that we only Legendre-transformed sources coupled to 0; i.e., the spinor 
field is treated as the environment. In Eq. ( p.28|) , the curved-spacetime Dirac 5 function 



is defined as in Comparison of Eq. ( |2.27| ) above with Eq. (4.13) of Ref. [|Ty] [which 



was computed for the 0{N) model] shows that the trlnFa^ in Eq. ( p.27| ) differs from the 
usual one-loop term by a factor of 2, owing to the difference (in the exponent) between the 
Gaussian integral formulas for real and complex fields |^ . 



The functional T2[(I),G, g^'^] can be evaluated in a loop expansion, which corresponds to 
an expansion in powers of h, 

oo 

T,[lG,gn = Y.h^r'^^\4>,G,gn, (2-36) 

1=2 

starting with the two-loop term, T^'^\ which has a diagrammatic expansion shown in Fig. |l|. 
The self-interaction leads to two terms in the two-loop part of the effective action, the 
second and third graphs of Fig. |l]. They are the "setting sun" diagram, which is O(A^), and 
the "double bubble," which is 0(A), respectively. The Yukawa interaction leads to only one 
diagram in F^^^ the first diagram of Fig. |l]. 



^aa'a" ^bb'b" 



d^x^/^ j d^x' ^/^'Gab{x,x')ii,^[Faib'{x,x')Fy'a"{x' ,x)], (2.37) 



where the trace is understood to be over the spinor indices which are not shown, and the 
three-index symbol c"'"^ is defined by 

{1 ii a = h = c = 
-1 if a = 6 = c = -, (2.38) 
otherwise. 

Here, we treat the A self-interaction using the time- dependent Hartree-Fock approximation 
pl| , which is equivalent to retaining the 0(A) (double-bubble) graph and dropping the O(A^) 
(setting sun) graph. We assume for the present study that the coupling A is sufficiently 
small that the O(A^) diagram is unimportant on the time scales of interest in the fermion 
production regime of the infiaton dynamics. The mean-field and gap equations including 
both the 0(A) and the 0{X'^) dia grams have been derived in a general curved spacetime in 
our previous paper |T^. 
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C. Evolution equations for <j) and G in curved spacetime 



The (bare) semiclassical field equations for tlie two-point function, mean field, and metric 
can be obtained from the CTP-2PI-CGEA by functional differentiation with respect Gab-, 
(pa, and g^'^, followed by identifications of and g'^'^ on the two time branches M, 



SS^ig^+mCg 



Sg'a' 

Sr[4>,G,g^^'^] 



6T[^,G,gn 



6G 



ab 



4>+—<t>- 



<t>+=<t'- 



<t>+—<i>- 



iiu nu ,,,, 

3+ =9- =9^ 



9^ =91 



(2.39a) 
(2.39b) 
(2.39c) 



Equations ( |2.39a| )-( |2.39c| ) constitute the semiclassical approximation to the full quantum 
dynamics for the system described by the classical action (|2.1|) . Equation (|2.39a]) should 
be understood as following after time branch indices have been reinstated on the metric 



tensor in the CTP-2PI-CGEA |10|. The field equation of semiclassical gravity (with bare 
parameters) is obtained directly from Eq. ( |2.39c| ). 



G,. + Acg,u + c^'^H^, + b'^^^H^, = -8nG{T^,), (2.40) 
in terms of the (unrenormalized) quantum energy-momentum tensor defined by 

2 6T[^,G,g>^'^] 



5g 



(2.41) 



7 7 7 M^^ n 



The energy-momentum tensor (T^^) is divergent in four spacetime dimensions, and must be 
regularized via a covariant procedure [^ ,[l|. 

Making the two- loop approximation to the CTP-2PI-CGEA, where we take r2 — t)?T'^'^\ 
and dropping the O(A^) diagram from the mean-field equation ( p.39b|) becomes 



A 



□ + + ^R{x) + -(p\x) + yG(x, x)j + hfii,^[Fab{x, x)] - h'g-^E{x) = 0, (2.42) 

where G{x,x) is the coincidence limit of Gab{x,x'), and in terms of a function S(?/) defined 
by 



s(y) 



d Xa/— ^ d X 



-g'^G++{x, x')tr,p [F++(x, y)F++{y, x')F++{x', x)] 

-G-+{x, a;')tr,p [F_+(x, y)F++{y, x')F+^{x', x)] 
-G+4x,x')tT^, [F++{x,y)F+4y,x')F.+ {x',x)] 

+G-- {x, a;')tr,p (x, y)F+_ {y, x')F_„ {x, x)] 



(2.43) 
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Making use of the curved spacetime definitions of the scalar and spinor field Hadamard 



kernels 36 



G^'\x,x') = {n\{^,ix),ipuix')}\n), 
F^'\x,x') = {n\[^M,^nix')m, 



and retarded propagators 



GR{x,x') = te{x,x'){n\ [ipn {x),ip^{x')]\n), 
Fnix^x') = te{x,x')m^M,^nix')}\n), 



(2.44a) 
(2.44b) 



(2.45a) 
(2.45b) 



the function S(?/) can be recast in a manifestly real and causal form, 
E(|/) = -2 j d^xy/^ j rfVv^Retr.p {e{x,x')G^^\x' ,x)F^^\x,x') 



-GRix,x'yFRix,x'))Fniy,x'yFniy,x) , (2.46) 



from which it is clear that the integrand vanishes whenever x or x' is to the future of y. The 
"gap" equation for Gab is obtained from Eq. ( p.39c| ), 



(G-i)'"(x, x') = ^'"(x, x') + ^c^"G(x, x)Six - x')^i= + hfc'"''''" €"''"11^, [F,,,,(x, x')F,»,.(x', x)] 

2 V-g' 

(2.47) 

Multiplying Eq. ( [^.47D through by Gab, performing a spacetime integration, and taking the 
++ component, we obtain 

'D + m' + ^R+ ^4>' + yG(x, x)) G++(x, x') 



+hf [ dx" J^'lC{x,x")G++{x\x') = -i5{x - x')^=, (2.48) 



(2.49) 



in terms of a kernel /C(x, x") defined by 

IC{x,x') = — itr^p [F++(x, x')F+4.(x', x) — F+_(x, x')F__|_(x', x)] . 
Making use of Eqs. ( |2.44b| ) and p.45b| ), this kernel takes the form 

/C(x, x') = Re tr,p [Fr{x, x')F^^^ (x', x)] , (2.50) 



which shows that the gap equation (|2.48| ) is manifestly real and causal. As will be shown 
below in Sec. |T| (in a perturbative limit), the kernel /C(x, x') is dissipative, and it reflects 
back reaction from fermionic particle production induced by the time-dependence of the 
infiaton variance. The gap equation ( |2.48| ) is therefore damped for modes above threshold^. 



^See P] for a similar discussion in the context of spinodal decomposition in quantum field theory 
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and this damping is not accounted for in the IPI treatments of inflaton dynamics (where 
only the inflaton mean field is dynamical). In contrast to previous studies [0,^,^ which 
assumed a local equation of motion for the inflaton propagator, the two-loop gap equation 
obtained from the CTP-2PI-CGEA includes a nonlocal kernel, which is a generic feature 
of back reaction from particle production. As stressed above, the dissipative dynamics of 
the inflaton two-point function can be important when the inflaton variance is on the order 
of the square of the inflaton mean-field amplitude; such conditions may exist at the end of 
preheating. 

The set of evolution equations ( ^2.42 ) for (p and ( 2.48| ) for G, is formally complete to two 
loops. Dissipation arises due to the coarse-graining of the spinor degrees of freedom. These 
dynamical equations are useful for general purposes, and are valid in a general background 
spacetime. However, in order to get explicit results, one needs to introduce approximations, 
as we now do. 



III. DYNAMICS OF SMALL- AMPLITUDE INFLATON OSCILLATIONS 

The effective evolution equations for the inflaton mean field and variance {(f"^) derived 
in the previous section are useful for studying fermion production when 00, the amplitude of 
the spatially homogeneous inflaton mean-field oscillations, is large, and the inflaton variance 
is of the same order-of-magnitude as (0o)^- As discussed in Sec. |I| above, such conditions 
can prevail at the end of the preheating period in chaotic inflation with unbroken symmetry 
Because of the dissipative kernel }C{x,x') in the gap equation ( p.48| ), which damps 
the evolution of G, and the back reaction terms in the mean-field equation, which damp the 
oscillations of 0, eventually the condition, 

/0o«m (3.1) 

will hold, at which point it is justifiable to follow the mean-field dynamics using the per- 
turbative, IPI, coarse-grained effective action [^]. Although in principle one should study 
this process in a general curved spacetime, for simplicity we assume spatial homogeneity, 
and that the inflaton mass is much greater than the Hubble constant, m ^ H. While this 
condition alone is in general not sufficient to ensure that curved spacetime effects are negli- 
gible during reheating (see, for example, 0, where cosmic expansion does affect preheating 
dynamics even though m ^ H), with the additional assumption of condition (^TTj) it is 
reasonable to neglect the effect of cosmic expansion in the spinor propagators p4| , pl| . In 
this and the following section, we also neglect the self-coupling A, because for the case of 
unbroken symmetry, the lowest-order A-dependent contribution to the perturbative inflaton 
self-energy is 0{fi^) [^, and we are only concerned with one- loop dynamics in this section. 



Let us therefore specialize to Minkowski space, and implement a perturbative expansion 
of the CTP effective action in powers of the mean field 0. This formally entails a solution of 
the gap equation ( |2.48| ) for G, a back-substitution of the solution into the CTP-2PI coarse- 
grained effective action, and a subsequent expansion of this expression (now a functional of 
only) in powers of 0. The resulting perturbative expansion for the effective action contains 
only free-field propagators. For consistency, one should use an initial density matrix for the 
spinor degrees of freedom which corresponds to the end-state particle occupation numbers 
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of the nonperturbative dynamics of Sec. |T[ For simplicity, however, we assume the initial 
quantum state for the spinor field is the vacuum state. Hereafter, Fab denotes the free- 
field, Minkowski-space, vacuum spinor CTP propagator, whose components are given by 
11,11311 



^,^0 = 7THie-^(--) , (3.2a) 



J {2tt)^ — fi^ + le 

F-4x,x') = - [ ^e-^(-') 5^ + "^^. , (3.2b) 
J [271)^ — fi^ — le 

F-+ix,x') = J ^e-^(---')2vr(i» + m)5(p2 -/i2)^(/), (3.2c) 

F^.{x,x') = I ^e-^(— ')27r(i, + m)5(p2 -/i2)e(V). (3.2d) 

The ++ and propagators admit a representation in terms of a time-ordering function 

eix,y) = eix^-y^), 

F^+{x,x')= 9{x,x')F^+{x,x') + 6'(x', x)F+_(x, x'), (3.3a) 
F_„(x,x')= e{x,x')F+^{x,x') + e{x',x)F.+ {x,x'). (3.3b) 

The CTP effective action can be expanded in powers of and we find 

r[0] = S'''[(i)] - ylndet(^'^^)^^ - thlndetFab + ri[0], (3.4) 

where the kernel A is defined by 

lA'^^x, x') = -c'^^D, + m^)5{x - x'), (3.5) 

and Fi is defined as —ifi times the sum of all one-particle-irreducible diagrams constructed 
with lines given by hA~^ and fiFah^ and vertices given by S^\^^'^ ^'^'\Jh and ?/','?/']/ fi. 

Because the free-field propagators A~^ and Fah do not depend on 0, the log(det) do not 
contribute to the variation of F[0] with respect to 0, and therefore, they can be dropped. 
The functional Fi[0] can be expanded in powers of /i, 

oo 

Fi[0] = ^;^'F«[0], (3.6) 
i=\ 

where the term F*^')[0] is the sum of all IPI /-loop graphs. Order by order in the loop 
expansion and the coupling constant, the CTP IPI effective action must satisfy the unitarity 
condition 

ril,,=,; ,_=, = o, (3.7) 

which has been verified to two-loop order in the case of scalar field theory |Q. The 
one-loop term in the loop expansion of the CTP effective action, F*^^''[0], can be further 
expanded in powers of 
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(3.8) 



n=l 



which corresponds to the usual amphtude expansion of the CTP effective action |2T1]. Fig- 
ure 1^ shows the diagrammatic expansion of F*^^). The terms F2^^[0] are generally divergent, 
but since the theory is renormalizable in the standard "in-out" formulation, it is renormal- 
izable in the closed-time-path, "in- in" formulation |]4^,|5|. 



A. One-loop perturbative effective action at 0(/^) 

The 0{p) term in the expansion of the one- loop CTP effective action, which is the first 
term in Fig. 0, takes the formQ 



^(l)r 



2 



d^xd'^x' (j)a{x)(j)a'{x')Fi,y{x, x')Fc'cix', x) . 



Making use of sum and difference variables 

1 



0+(x) + 0_(x') 



the functional Fn^-*! 



A{x) = (f)+{x) — 0_(x' 
can be recast in the form 



(3.9) 



(3.10a) 
(3.10b) 



d'^xd'^x' 



T.{x)A{x')V2{x,x') + -A{x)A{x')U2{x,x') 



in terms of manifestly real kernels 1^2(2;, x') and ^2(0;, x') defined by 

V2{x,x') = Imtr,jF++(a;,x')F++(a;',x) + F+_(x,a;')i^-+(a;',a;)], 
Af2{x,x') = — Retrsp[F++(x, x')F++(x', x)]. 



(3.11) 



(3.12a) 
(3.12b) 



Only the kernel V2{x,x') contributes to the mean-field equation of motion. The kernel 
J\f2{x,x') constitutes a correlator for noise, and will be discussed in Sec. |V[ The unitarity 
condition ( |3.7| ) requires that the sum of diagrams proportional to S(x)S(x') vanish identi- 
cally. With the definitions of the retarded spinor propagator, Eq. (p.44b| ), and the spinor 
Hadamard kernel, Eq. ( |2.45b| ), which in Minkowski space take the form 



Fji{x,x') = i9{x,x')[F^^{x,x') — F+_(x, x')], 
F^^) (x, x') = F_+ (x, x') + F+_ (x, x') , 

the kernel T>2{x, x') can be written in a manifestly causal form. 



(3.13a) 
(3.13b) 



^Note that there are no nonzero graphs with an odd number of vertices in this model. 
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V2{x,x') = ]^Reii^^[Fn{x,x')F^^\x\x)]. (3.14) 

Using Eq. (|3.14|) , one can compute V2{x,x') in an arbitrary curved background spacetime. 
It should be noted that T>2{x^x') is just the lowest-order term in the series expansion of 
K.{x,x') [defined in Eq. ( p.50| )] in powers of the coupling constant /. The appearance of the 
retarded propagator in Eq. ( |3.14| ) guarantees that the contribution of r^^^ to the mean-field 



equation of motion is causal. 

Let us now evaluate V2{x,x') using dimensional regularization and the modified mini- 



mal subtraction (MS) renormalization prescription |^5|. Dimensional regularization requires 

changing the coupling constant so that the interaction has the correct dimensions in n 
spacetime dimensions, 

/ /A(^-")/^ (3.15) 

where we have introduced a parameter A, the renormalization scale, which has dimensions 
of mass. By Lorentz invariance and causality, the product of Feynman propagators can be 
written in terms of an amplitude A2 [^], 

tT^AF++{x,x')F++{x',x)]=^ I ^^e-^^^--'^A2{e + ie), (3.16) 

and with this choice of renormalization prescription, the amplitude A2{k'^) takes the form 

A2{k^) = daE{a; P)log (^^^) , (3.17) 

where E{a] k"^) is defined byQ 

E(a; A;2) = _ a(i _ a)fc2. (3.I8) 

Note that in Eq. ( |3.17| ), the a integration shows up via the Feynman identity 
1 

C,...C 



N 



(A^— 1)! / dai--- I daN5{ai + ■ ■ ■ + aN-i)[oiiCi + ■ ■ ■ + unCn] ^- (3.19) 



The ie appearing in Eq. ( |3.16| ) ensures that the amplitude A2 is evaluated on the physical 
sheet . The logarithm in Eq. ( p.l7| ) has a negative real argument when the two condi- 

tions /c^ > 4/1^ and |2a— 1| < A/r^-4/?7p are both satisfied. When |2a;— 1| < a/1 — ifi'^/k'^, 
the amplitude A2{k'^) has a branch cut (considered as an analytically continued function of 
k^) for (/c°)2 > k'^ + 4/i2. The discontinuity across the branch cut is related to the "cut" 
version of the diagram [the second term in Eq. (|3.12a|) ] via the Cutosky rules [^ , ^ ,^8|-|5T 



^The notation E{a; k"^) used here should not be confused with E{k), the complete elliptic integral 
of second kind. 
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tTjF+4x,x')F_+{x',x)] = -i j 



(2^ 



-ik{x—x') 



Disc[A2(fc2)]^(A;0). 



(3.20) 



From Eqs. ( |3.16D , ( |3.17| ), and ( |3.20[ ), it is straightforward to obtain an expression for the 
dissipation kernel, 



V2{X,X') 



47r2 J (27r 



log 



\E{a;k^) 



in9[-E{a; k'^)]sgn{k°) 
(3.21) 



where we have now taken the limit e ^ 0+. One can verify by inspection that this kernel 
is real. However, the second term in Eq. ( p.21|) breaks time-reversal invariance and leads to 
dissipative mean field dynamics. The one-loop Fourier-transformed mean-field equation is 
(dropping the caret from 0) 



k-m^ + ik"^2{k) 



3^ 

47r2 



daE[a; k jlog 







A2 



0(A;) = -J{k), (3.22) 



where 72(fc) is the dissipation function, defined as —ihp/k^ times the Fourier transform of 
the second term in Eq. (|3.21J ), 



Stt \k^\ 



k^ 



(3.23) 



The one-loop 0{p) dissipation kernel agrees with previous calculations of the probability 
to produce a fermion particle pair of momentum k In Eq. ( |3.22| ), J{k) is an 



external c-number source. The imaginary term ik^'^2'P in Eq- (|3.22|) breaks time-reversal 
invariance and acts as a A;-dependent dissipative force in the mean field equation. The 9 
function enforces the energy threshold for the virtual fermion pair in the one- loop 0{p) 
diagram to go on-shell. The dissipative mean-field equation ( p.22| ) is essentially the linear- 
response approximation to the effective inflaton dynamics. It should be noted that the 
dissipation kernel V2{x,x') is nonlocal, in contrast with the local friction-type dissipation 
assumed in earlier studies of post-inflation reheating [Q. However, in the limit fi^ 0, 
the dissipation kernel does become time-local, as there is no longer a length scale in the 
expression for V2{x,x') which could define a time scale for nonlocal dissipation M. 



B. One-loop perturbative effective action at 0(/^) 

The 0{f^) term in the one-loop OTP effective action consists of the "square" diagram, 
which is the second term in Fig. ^ 



^(l)r 



,abc^a'b'c'^def^d'e'f' / d\d' x' d'y' d'y ir^,[F,y {x, x')F, f, {x\ y')F,, e{y\ y)Fj,{y , x)] 



X^a{x)^a'ix')(f)d{y)(i)d'iy') 



(3.24) 
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Expanding out the contracted CTP indices, we obtain 



^(l)r 



d xd X d y d y 



I 
4 

0+(x)0+(x')0+(y)0+(y')trsp{i^++(x, y')F++{y\ y)F++{y, x)} 

+0_(x)0_(xO0-(y)0-(l/Otrep{F--(x,x')F__(x',i/')i'--(l/',Z/)i'-~(z/,a;)} 
-40+(x)0_(xO0-(yO0-(z/)trsp{F+-(x,x')F__(x',y')^--(z/',l/)^-+(l/,a;)} 
-40_(a;)0+(a;O0+(yO0+(y)trsp{F-+(x,x')F++(x',y')^++(y',y)^+-(2/,a;)} 

+40+(x)0+(xO0-(yO0-(l/)tr.p{F++(x,x')F+^(x',i/')^--(l/',Z/)^-+(y,^)} 

+20+(x)0_(xO0+(y')0-(i/)tr.p{F+-(x,a;')F_+(x',y')^+-(2/',y)^-+(y,^)} 



(3.25) 



When r4^'' is expressed in terms of A and E [defined in Eqs. ( p.lOa|) and (|3.10b|) ], only those 
terms with one factor of A and three factors of S contribute to the mean field equation of 
motion. As a consequence of the unitarity condition (|3.7|) , the sum of terms proportional 
to four factors of S must vanish. Keeping only those terms in the effective action which 
contribute to the mean field equation or are quadratic in A, we find 



d^xd^x'dyd^y A{x)J:{x')J:{y')^{y)V4{x, x', y', y) 



+ -A(x) A(x')S(y')S(2/)M(x, x\ y\ y) 



(3.26) 



in terms of a kernel I^4(x, x', y', y) defined by 
Vi{x,x',y\y) = -Imtr 





[x, X 


')F++ 


{x',y')F^+ 


{y'.y)F++ 


iy,x) 


+F^ 


_i_(x 


x')F+ 


-{x\y')F_ 


-{y\y)F- 


+iy,x] 


+F-, 




x')F. 


-ix',y')F^ 


+{y',y)F, 


-+{y^ Xj 


+F-, 




x')F^ 


+{x\y')F^ 


Ay',y)F^ 


+iy,x] 


-F^ 




x')F_ 


+ix',y')F, 


+{y',y)F, 


-+i.y,x] 


-F^ 


+ {x 


x')F+ 


~-{x',y')F- 


+{y\y)F^ 


-+{y,x] 


-F^ 


+ [x 


x')F+ 


+ ix',y')F, 


-{y',y)F- 


-+i.y,x] 




_(x 


x')F. 


-ix',y')F- 


-{y\y)F- 


-+{y,x] 



(3.27) 



and a "noise" kernel A/'4(x, x', y', y) defined by 
X4:{x,x',y',y) = Retr 



F+4x 


,x')F++ 


{x',y')F++ 


iy',y)F++iy,x) 


+F++ 


[x, x')F^ 


-{x',y')F_ 


-{y',y)F^+{y,x) 




[x, x')F^ 


-{x',y')F- 


-+{y\y)F++{y,x) 




[x, x')F_ 


+ix',y')F, 


-^{y\y)F^+{y,x) 


+F+^ 


[x, x')F_ 


+{x',y')F^ 


^+{y',y)F++{y,x) 


-F++ 


[x, x')F+ 


-{x',y')F_ 


.+ iy',y)F++{y,x) 
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x')F++ 




{y',y)F^+{y,x) 










iy',y)F-^ 


Xy,x) 




Ax,y')F+ 


4y',x')F+_ 


,ix',y)F^ 


+{y^x) 






y')F- 


{y',x')F^+ 


{x',y)F+^ 


-iy,x) 






y')F+- 


{y',x')F__ 




-iy,x) 


+F+- 


(a;, 


y')F^+ 




(x',?/)F_4 


Xy,x) 


-F+- 




y')F-+ 


{y\x')F++ 


{x',y)F+^ 


-{y,x) 






y')F+- 


{y\x')F_+ 


{x',y)F+^ 


Xy^x) 




(x, 


y')F++ 


iy\x')F+^ 




-iy,x) 






y')F- 


{y',x')F__ 


{x',y)F_^ 


Xy,x)} 



}/2 



(3.28) 



The noise kernel ACt does not contribute to the mean field equation of motion. There are, 
of course, terms in r4^''[0] which are higher order in A, for example, O(A^), but in passing 
over to a stochastic equation for (p in Sec. [TV] , we will be assuming that A is small, so that 
higher-order terms in powers of A can be ignored. Such terms will in general contribute to 
non-Gaussian noise, which will be studied in an upcoming paper 

Let us evaluate the first term of Eq. ( p.27|) , which consists of only Feynman propagators. 
The term is logarithmically divergent, and as in Sec. [ill A| , we use dimensional continuation 
and the modified minimal subtraction (MS) renormalization scheme. Because we are only 
interested in deriving the dissipative terms in the mean-field equation coming from this 
diagram, and because we are assuming m ^ /i, we include only the one-loop logarithm. We 
find 



tr. 



F++{x, x')F++(x', y')F++{y', y)F++{y, x) 



d*h d^k2 d% 



(27r)4 (27r)4 (27r)4 
where the amplitude A4{ki, k2, k^) is defined by 



- log only 

g-ih(fel+fc2+fc3)-X+fclV+fc2V+fc3-?/]^^(^^^ ^^^^ (3^29) 



A^iki, k2, k-i 



27r2 



daida2da^ log 



-E4(ai,a2,Ct3;fci,fc2,fc3) 
A2 



(3.30) 



in terms of a function defined by 

Ei{ai,a2,a^]ki,k2,k^) = 



1 - ai)ki (1 - «! - a2)k2 + {1 - ai - a2 - 0(3)^3 

-(1 - ai)kl - (1 - ai - a2)ikl + 2ki ■ /ca) 

-(1 -ai-a2- a3){2ki ■ ^3 + 2k2 ■ k3 + k^) + (3.31) 



As in Sec. [Ill A] , the cut diagrams in Eq. (|3.27|) are related to the discontinuities in 
i?4(ai, a2, «3; ^1, ^2, k3) via the Cutosky rules. The details are shown in the appendix. We 
can then express the dissipation kernel V4{x,x',y',y) as a Fourier transform over external 
momenta, 
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d^ki d^k2 d^k^ 



(27r)4 (27r)4 (27r)4 
in terms of a function 'D4^{ki, k2, k^) defined by 

3 



'D4(/ci, k2, ks 



2n 



1 /" , , , , f\E4^{ai,a2,as;ki,k2,ks) 
- \ daida2da3 log I — 

+^sgn(A;° + k^)h[{k2 + k^f] + zsgn(A;? + k^^)h[{k, + fcs 
+isgn(A;° + k°2 + ^3)^[^2] + isgn(fc°)/i[(A;i + k2 + k^f] 

+isgn{k^)h[kl] + isgn{k^^)h[kl] - iH{ki, k2, k^ 



(3.33) 



and the functions h{s) and H{ki, k2, ks) are defined by 



H{ki, /C2, /cs) = / {9[-Ei{ai, 0^2, 03; ki, k2, k^)] 

'ai,O2,O3>0 



(3.34a) 



(3.34b) 



Equation ( 3.33 ) leads to the following mean-field equation at 0(/'^), 



k' - m' + ik"^2ik) - 



47r2 



daE{a; k )log 



A2 



0(fc) 



3i^ 
27r 



(27r)4 (27r) 



— f daida2da3\og ( - — ~^ ^ — — li 

TT J \ 

+sgn(g° + f)h[iq + l)^] + sgn(A;°)/i(g2) 

+sgn(A;° - f)h[{k - If] + sgn(g°)/i(A;2) 

+sgn(A;° - g° - f)h{P) 

+sgn{l')h[{k-q-lf] 



+H{l + q- k,-q, -I) 



-J{k). (3.35) 



The presence of terms of the form isgn{p^)h{p'^) in Eq. (|3.35| ) clearly signifies dissipative 
dynamics. The ^-function in Eq. ( ^.34aj ) enforces the energy threshold for the virtual fermion 
quanta created at a particular vertex to go on-shell. Comparing Eq. ( |3.35| ) and Eq. ( p. 221) , 
and assuming spatial homogeneity, we see that the 0{f^) dissipation kernel must be taken 
into account whenever the condition ( |3.1| ) fails to hold for the solution 0(t) to Eq. (p.22| ). 

At the end of the regime of parametric resonance in chaotic inflaton, i.e. the "preheating" 
regime, the inflaton mean field may oscillate with an amplitude as large as ~ m/(yf^^, where 
g^-^ is the coupling to another scalar field x, typically on the order of 10~^ @]. Condition ( p.l|) 
would then be violated if / > g^^. In this case it would be necessary, at a minimum, to take 
into account higher order terms (such as P4) in the mean-field equation, until such time as 
the amplitude (po{t) has decreased to the point where Eq. is satisfied. 
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IV. NOISE KERNEL AND STOCHASTIC INFLATON DYNAMICS 



Although the kernels Af2{x,x') and N'i{x,x' ,y' ,y) do not contribute to the mean field 
equation, i.e., the equation of motion for 0, they contain information about stochas- 
ticity in a quasi-classical description of the effective dynamics of the inflaton field 
|, |31|j30| , |32|j55| -p^,p!^ . In this section, we study the effect of stochasticity on the dynamics 



of the inflaton mean fleld, within the perturbative framework established above. 



A. Langevin equation and fluctuation-dissipation relation at 0{f^ 



In this section we show how to obtain a classical stochastic equation for the inflaton 
fleld from the 0{f^) perturbative CTP effective action. From Eq. ( p. Ill ), it follows that the 
0{f^) one-loop perturbative CTP effective action has the form 



T{(f)] =S^\(t)\ + / (fxd^x' 



S(x)A(x')/i2(a;,x') + -/\{x)A{x')p2{.x,x') 



where for simplicity we have deflned 



1'2{X,X') = hf'^J\f2{x,x'), 
/i2(x, x') = hf^V2{x,x'). 



(4.1) 



(4.2a) 
(4.2b) 



In order to extract the stochastic noise arising from the kernel A/2(x, x'), we use the Gaussian 
identity 



exp 



-— / d'^xd'^x'A{x)A{x')iy2{x,x') 



I 



N / D^exp 



— J d^xd^x' i2{x)i'2 ^{x-,x')i2{x') 
+ \ I d'xUx)A{x)\, (4.3) 



where is a normalization factor which does not depend on A, and ^2 is a c-number 
functional integration variable. Following [^], we now deflne a functional 



P[6] = AT exp 



h 



d^xd^x' i2{x)v2^{x,x')i2{x') 



(4.4) 



and it follows from Eq. ( |4.3| ) that P[^2\ is normalized in the sense of 

D^2PU = 1. (4.5) 

Using Eq. ([4.3|), we can rewrite the 0{p) one- loop CTP effective action, Eq. ( |4.1| ), as 

T[4)] = -zh\og j DG^[G]exp + j d'^xd^x'^{x)A{x')fi2{x,x') + j d'^x^2{x)A{x) 

(4.6) 
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This suggests defining a new effective action which depends on both ^2 and 0± (dropping 
the carat from 6), 



Let us define a type of ensemble average 



and note that Eqs. 



and ( |4.8| ) imply that 

((6(x))) = o, 

{Mx)Ux'))) = hu2{x,x'). 



Clearly then, 



-{{mm 



't>+—4>- 



-m 



<t>+—4>- 



Taking the variation of r[(/), ^2] with respect to 0+ and setting 
(after a Fourier transform) 



(4.7) 



(4.8) 



(4.9a) 
(4.9b) 



(4.10) 



i), we obtain 



47r2 



2 /•! 



daE{a; k^)\og 



\E{a;k'^ 
A2~ 



4>(k)+Uk) = -Jik), (4.11) 



where ^2(^) is defined by 



Uk) = I d^xe'^'^Ux). 



(4.12) 



We now interpret Eq. ([4.11| ) as a Langevin equation with stochastic force ^2- The infiaton 
Fourier mode appearing in Eq. ( 4.11 ) should be viewed as a c-number stochastic variable, 
and the presence of the stochastic force ^2 indicates spontaneous breaking of spatial transla- 
tion invariance by a Gaussian (but not white) noise source .^2 1211 • Moreover, this stochastic 
equation obeys a zero-temperature fluctuation-dissipation relation, as we now show. First, 
let us calculate the one-loop 0{p) noise kernel, Af2{x,x') [defined in Eq. ( p.l2b|) above], 
using dimensional regularization and modified minimal subtraction. 



1^2(3;, X 



in 



d^k 



-ik{x—x') ~ 



i^2ik), 



in terms of the Fourier-transformed noise kernel 



Stt 



fc2 



(4.13) 



(4.14) 



The noise kernel i'2{x,x') is colored; colored noise has been observed in other interacting 
field theories pT|j3^ , |55[] . By inspection of Eqs. ( |3.23| ) and ( [4.14| ), it follows that 
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|fcO|72(A;) = Uk), (4.15) 
which leads to the zero-temperature fluctuation-dissipation relation 



U2{t,k)= dt'K{t~t')-f2it',k), (4.16) 



in terms of the distribution-valued kernel K{t) defined by 

/•OO 7 

K{t) = / — cucos(c^t), (4.17) 

Jo TT 

and the spatially Fourier-transformed dissipation function and noise kernel, 

Mt,k)= / —e''''"'D2{k), (4.18a) 



-oo 27r 

72{t,k)= I -^e-^''%{k). (4.18b) 

This shows the physical significance of the noise kernel i'2{x,x') in an effective description 
of the dynamics of the scalar mean field. 



B. Langevin equation and fluctuation-dissipation relation at 0{f'^) 

In this section we consider the 0{f'^) one-loop noise kernel, A/4. The non-normal- 
threshold singularities in lead to a noise kernel which depends on S, which is known 
to lead to ambiguities in the resulting Langevin equation [^,^. The meaning and inter- 
pretation of the non-normal-threshold parts of A/4 and P4 will be the subject of a future 
study p3. Here, we focus on the effect of the normal-threshold singularities of ^4, which 



for the noise kernel, A/4, contribute a term 



-j rf^xrfVA(x)A(x')S(x)E(x)z/4(a;,x') (4.19) 
to the CTP effective action, where the kernel 1/4(0:, x') is defined by 

^4(x,x') = / 7?!ie"^'-^^"^'^M9'), (4.20) 



vr J (27r)4' 

and the function h{s) was defined in Eq. ( |3.34a| ) above. The normal-threshold singularities 



of the dissipation kernel, 1^4, [the second and third terms of Eq. ( |3.33| )], lead to the following 
contribution to the CTP effective action, 

j ciWa;'A(x)S(x) [S(a;')]V4(a;,x'), (4.21) 

where the kernel ^^{x^x') is defined by 
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H4{x, x') 

With the definitions 

/i4(a;, X' ] = I 
z/4(x, X 



3ihf^ f dS 



-iq {x'-x) ( 



vr J (27r)^ 



sgn{q )h{q 



-iq-{x' —x) ~ 



(4.22) 

(4.23a) 
(4.23b) 



it follows immediately that the normal-threshold parts of and A/4 obey a fiuctuation- 
dissipation relation identical in form to Eq. (|4.15|) , 



Q 174(g) = I'iiq)- 



(4.24) 



Making use of Eq. ( [4.3|) , the 0{f^) effective action (including only normal-threshold contri- 
butions) can be written in the form 

V[(P,i2.iA\=S^[(P] + j dWx'A(a;)S(x')/i2(a;,x') + j dSUx)A{x) 

+ j d^xd^x'A{x)J:{x)[J:{x')f fi4{x,x') + j d^x^i{x)A{x)J:{x), (4.25) 



where the stochastic noise source ^4 satisfies the conditions 

{{Ux)Ux'))) = hiy,ix,x'). 



(4.26a) 
(4.26b) 



Taking the functional derivative of Eq. ([4.25|) and making the usual identification, we obtain 
a Langevin equation with an additive noise ^2 and a multiplicative noise ^4, 



q - m +iq 72(g) - 



4772 



2 /•! 



daE{a] q^) log 







\E{a-.q')\ 



27r2 



d^a log 



\Ei{ai, ^2, as; / - g, -k, I - k)\ 



A2 



-6(g) -J(g) 



(2vr) 



(4. 



where d^a = daida2da3. The stochastic force 1^4 is clearly seen to contribute multiplicatively 
to the Langevin equation for (p. 



C. Homogeneous mean field dynamics at 0{f'^) 

To make connection with post-infiationary reheating, it is customary to assume that the 
mean field is spatially homogeneous In this case, the Langevin equation ([4.1 1|) 

takes the form 
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[u^ -m^ + iul3{u) + r/(cu)] ^{uj) + |2(^) = -J(c^), (4.28) 
where we have defined 

riiu^) = -^-K f daE{a; u') log ( \.E1^2^\ . (4.29b) 



The total energy dissipated to the fermion field over the history of the dynamical evolution 
of the mean field is given [at 0{py\ by 

S = -J^dtF^{t)^, (4.30) 

where the friction force F^{t) is the Fourier transform of if3{uj)uj(f){uj). After a bit of Fourier 
algebra, we obtain an expression for the ensemble-averaged, total dissipated energy. 



It is straightforward to compute the variance in the total dissipated energy. Making use of 
Eq. ( Pb| ), we find 



(4.32) 



where the function /(cu) is defined by 



/.^c^2_4^2 . 2 \3/2 

/(u;) = y dkeicu' - e) M - -rz:^] ■ (4-33) 

Following |T^, we assume that the infiaton field is held fixed via an external, constant 
c-number source J for t < 0, and that the source is removed for t > 0, 

J(t) = je{-t). (4.34) 

Setting h = 1, assuming that m S> /x, and expanding to lowest order in /, we obtain for the 
ensemble averaged dissipated energy, 

m = -r^r (4.35) 

Let us now compute the variance in the total dissipated energy. Performing a regularization 
via dimensional continuation, we obtain 
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(4.36) 



where 5 is a constant of order unity defined by 5^ = |119/60 — 7em — log(47rm^/A^)|, and 7em 
is the Euler-Mascheroni constant, ^ 0.5772. Taking the ratio of the square root of Eq. ( [4.36| ) 
and Eq. ( |4.35| ), we obtain the relative strength of the rms fiuctuations in the total dissipated 
energy density, £,^,, 

m) m) v^vrj- ^ ■ ^ 

The parameter J is related to the initial inflaton amplitude 0o(io) by J = 0o(^o)'^^/2, which 
leads to 

'^"^^ - 0.390^. (4.38) 



{{£)) v^7r0o(to) ' Mto) 

The fundamental assumption which justified the perturbative expansion in /, Eq. ( |3.1| ), is 
seen to be independent of Eq. (|4.38| ). Therefore, the ratio £,ms/£ is not required to be small 
by consistency with perturbation theory. As the initial infiaton amplitude 0o is made larger, 
the relative strength of the rms fiuctuations of £ is seen to decrease, in accordance with the 
correspondence principle. It has been shown that the fiuctuations in the total dissipated 
energy density are related to the fiuctuations in the occupation numbers of modes . 

Let us now examine whether the rms fiuctuations in the total dissipated energy, as given 
by Eq. ( [4.38| ), is significant, given a reasonable value for the infiaton amplitude at the end of 
the preheating regime (the period of parametric resonance- induced particle production). In 
chaotic infiaton with a scalar field x coupled to the infiaton field via a coupling constant g^^, 
the typical infiaton amplitude at the end of the preheating regime is on the order of m/g^^ 
In this case, we would find £rmj {{£)) — fg<j>xy from which it is clear that fiuctuations 
in the total dissipated energy are not significant relative to the ensemble-averaged total 
dissipated energy, and therefore should not appreciably affect the reheating temperature. 
However, in new inflation scenarios where the infiaton amplitude 0o can be on the order of 
m at the onset of reheating, the ratio m/0o can be of order unity ^3|- ^^is case, the ratio 
£vn^J {{^)) — f 1 which may not be a negligible effect. 

Although as shown above, stochasticity does not dramatically affect the total energy 
dissipated via fermion production in chaotic infiation, we may inquire whether the noise term 
in the Langevin equation for the infiaton zero- mode, Eq. ( [4.28|) , may nonetheless be non- 



negligible during the reheating period. Let us compute the rms fiuctuations in the infiaton 
zero-mode, 0rms- Following methods described above, we find that the rms fiuctuations of 
the infiaton zero mode are given, to 0(/^), by 

-M, (4.39) 



TT V60 



TT 



where cr^ = |61/30 — 7em — log(47rm^/A^)|. Equation ( [4.39|) is seen to be independent of the 
infiaton zero- mode amplitude 0o- 



28 



In order to determine the relative importance of fluctuations in the inflaton zero-mode 
amplitude 0o during and at the end of the reheating period, we must introduce curved 
spacetime arguments. This is because the end of the reheating period is determined by 
the time t^.^^ at which the Hubble constant becomes of the order of 3/3(m) for the case 



of A = being discussed in this section ||53|,|23[|. Starting with the semiclassical Einstein 



equation ( p. 401) for spatially flat Friedmann-Robertson- Walker (FRW) cosmology, setting 
b = c = Ac = (following arguments similar to those of Sec. Ill D of Ref. |]I|]), retaining only 
the inflaton zero mode as the dynamical degree of freedom (for consistency with FRW), and 
retaining both the ip field energy density and the classical, stochastic energy density of 
the inflaton zero mode, we have 



where a is the scale factor, the dot denotes a derivative with respect to cosmic time, and 
p(02) is the energy density as a function of the time-average (over one period of oscillation) 
of 0^, which is given by the virial theorem, 

p{4>) ~ m202 = lrn2(0o)2. (4.41) 



Making use of Eqs. ( [4.28|) , (|4.34|) , and (|4.40| ), we obtain an approximate expression for the 



(ensemble-averaged) energy density of the inflaton zero-mode at the end of the reheating 
period[] (to lowest order in /), 

P ^ \^ ; , 4.42 

(87r)'^e 

where e is the base of the natural logarithm. Note that this expression is independent of 
the initial inflaton amplitude 0]. Equation ( [4.42| ) allows us to solve for the value of (po at 
the end of reheating. We find 



Mt...) - ^(3^)3// • (4-43) 

The rms fluctuations in the inflaton zero mode, can only play a role in the inflaton 
zero-mode dynamics during reheating if the ratio 0,.ms/0o is not small relative to higher-order 
[e.g., 0(/^)] processes which we are neglecting. In light of the minimum inflaton zero-mode 
amplitude attained during reheating, Eq. ( [4.43| ), we find that the ratio of fluctuations in the 
inflaton zero-mode to the zero-mode amplitude is given by 



^ We wish to emphasize, however, that this expression does not take into account the regime 
of nonperturbative dynamics discussed in Sec. ||, and therefore should not be expected to yield 
a correct reheating temperature in a realistic inflationary scenario. However, it suffices for the 
present discussion of rms fluctuations of the inflaton amplitude, where we assume an idealized case 
similar to Eq. (106) of Ref. p|. 
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Sam je m 
7MpV45 "^-^ /Mp 



(4.44) 



We estimate the ratio of the mean-squared inflaton amphtude fluctuations to the shift in 
the inflaton mass to be 



()2 

rms 



ri{m) 



0.01. 



(4.45) 



If, prior to the end of reheating at t^^^, 0rms/0o(^) becomes larger than higher-order terms 
which are neglected in our perturbative expansion, then fluctuations in the inflaton zero 
mode are a non-negligible effect. This will happen when 0rms/0o(^ond), given by Eq. ( [4.44p , 
is not ^ 1. This shows that stochasticity must be taken into account in the dynamics of 
the inflaton zero mode, during the late stages of the reheating period. 



V. CONCLUSIONS 

In this paper, we present the results of a study of (unbroken-symmetry) inflaton dynamics 
during the late stages of reheating, which is dominated by fermion particle production to a 
light spinor field coupled to the inflaton field via a Yukawa coupling. We derived coupled 
nonperturbative equations for the inflaton mean field and two-point function, in a general 
curved spacetime, and showed that, in addition to the dissipative mean-field equation, the 
gap equation for the two-point function is also dissipative, due to fermion particle production. 
Simultaneous evolution of the inflaton mean-field and two-point function is necessary for 
correctly following the inflaton dynamics after the end of the preheating period, because the 
large value of the variance invalidates use of the ordinary perturbative, IPI effective action. 

We also derived the dissipation and noise kernels for the small-amplitude dynamics of the 
inflaton field, valid in the late stages of reheating when the inflaton mean-field amplitude is 
very small. The 0{p) noise and dissipation kernels, as well as the normal-threshold parts of 
the 0{f'^) noise and dissipation kernels, are shown to obey a zero-temperature fluctuation- 
dissipation relation. With the noise and dissipation kernels, a Langevin equation for the 
inflaton zero mode is derived, and it is shown that the noise leads to a variance for the 
inflaton amplitude which is non-negligible before the end of reheating. 
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APPENDIX A: DISCONTINUITIES OF THE SQUARE DIAGRAM 



In this appendix, the seven terms of Eq. ( ^.271 ) involving cut propagators are exphcitly 
evaluated using the Cutosky rules. The second and third terms of Eq. (|3.27|) correspond to 
normal-threshold singularities in the t and s channels, and are given by 



(2vr) 



F++(g)F+_(g + h)F_.{q + h + k2)F_+{q + k, + k2 + h 

I 



iDisc[Ai{ki, k2, /i;3)|ai=a3=o]6'(fc2 + ^3) 



(Al) 



and 



F+4q)F_4q + k,)F^+{q + k^ + k2)F++{q + k, + k2 + k^) ^ 

= -iDisc[A4{ki, k2, A;3)|a2=0; ai+a3=l]^(^l + ^2)5 



(A2) 



respectively. The third term in Eq. (p.27|) corresponds to the leading-order singularity of 



the square diagram (i.e., the solution of the Landau equations in which ai, 02,0(3 are 
all nonzero). 



d^q 
(2vr) 



+ ki)F+^{q + ki + k2)F^+{q + fci + ^2 + ^3) 

= zDisc[A4(A;i, k2, k3)\c.,^^,,c.,>o]0{kl)e{-kl)9{kl). (A3) 



The last four terms in Eq. ( p.27| ) correspond to the four remaining twice-contracted singu- 
larities, and are given by 



dS 



(2vr) 



4trap 



F++{q)F+_{q + ki)F_+{q + k^ + k2)F++{q + k, + k2 + k^ 



zDisc[yl4(/ci, ^2, k. 



3)\a2=as=oW{k2), 



(A4) 



d^q 
(27r) 



F++{q)F++{q + ki)F+^{q + fci + k2)F^+{q + ki + k2 + k^) 
= iDisc[A4(A;i, A;2, k^)\c,.,=o. „,+c.2=i]^(^3 



(A5) 



d^q 



(2vr) 



^tr^p 



F+_{q)F_+{q + ki)F++{q + ki + k2)F++{q + ki + k2 + k^) 
= zDisc[A4(A;i, ^2, k3)ia^=a2=oW{ki), 



(A6) 



d^q 



(2vr) 



4^^sp 



F+_{q)F.^{q + fci)F__(g + fci + k2)F^+{q + ki + k2 + k^) 

= zDisc[A4(A;i, k2, A;3)|ai=0; a2+as=i]0{ki + k^2 + ks). (A7) 
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FIGURES 




FIG. 1. Diagrammatic expansion for F*-^), the two-loop part of the CTP-2PI-CGEA. Sohd Unes 
represent the spinor propagator F (as defined in Sec. U), and dotted Unes represent the scalar 
propagator G. The vertices terminating three cj) lines are proportional to the scalar mean field (j). 
Each vertex carries spacetime (x) and CTP (+, — ) labels. 




FIG. 2. Diagrammatic expansion for F'-^^ the one-loop part of the CTP-IPI-CGEA. Solid lines 
represent the spinor propagator F (as defined in Sec. ID), and dotted lines represent multiplication 
by the scalar mean field cf). Each vertex carries spacetime (x) and CTP (+, — ) labels. 
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